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Abstract 

The statistical properties of fully developed hydrodynamic turbu- 
lence can be successfully described using methods from nonextensive 
statistical mechanics. The predicted probability densities and scal- 
ing exponents precisely coincide with what is measured in various 
turbulence experiments. As a dynamical basis for nonextensive be- 
haviour we consider nonlinear Langevin equations with fluctuating 
friction forces, where Tsallis statistics can be proved rigorously. 



1 Introduction 



The formalism of nonextensive statistical mechanics was introduced by C. Tsal- 
lis in 1988 and has been further developed by many others [0|-||6[. Recently, 
it has been pointed out in a number of papers that the nonextensive approach 
is in particular useful for the description of stochastic properties of fully devel- 
oped turbulent flows |]^]~ []l5|| . In this paper we outline the main idea underlying 
this new statistical mechanics approach to turbulence und compare with various 
experimental measurements. 

First we will summarize the basic idea of the nonextensive formalism and 
point out some interesting connection to systems with fluctuating temperature or 
fluctuating energy dissipation (i.e. nonequilibrium systems with a stationary state, 
see e.g. As an example where everything can be calculated analytically, 

we will consider suitable generalizations of Langevin equations with fluctuating 
friction forces in section 3 |17|. These underly the generalized statistical mechanics 
in an analogous way as ordinary Langevin equations underly ordinary statistical 
mechanics. The concrete application to turbulent flows will be described in section 
4. Finally, in section 5 we compare with various experimental data in turbulent 
flows. The agreement turns out to be very good, thus indicating that methods from 
generalized statistical mechanics are a useful tool for modeling turbulence. Fully 
developed turbulent flows appear to effectively extremize the Tsallis entropies. 



2 Tsallis statistics and fluctuations 

In the nonextensive statistical mechanics approach, the usual Shannon entropy 
Si = — X^jPilnpi is replaced by the more general Tsallis entropies 

The Sq are related but different from the Renyi entropies (see, e.g., The 
Pi are probabilities associated with the microstates i of the physical system, and 
g 7^ 1 is a parameter, the 'entropic index'. The ordinary Shannon entropy Si is 
recovered for g — > 1. 

Extremizing Sq under suitable constraints (see e.g. |^] for a review), just in an 
analogous way as it is done in ordinary statistical mechanics, one formally obtains 
a g-generalized version of the canonical distributions. The probabilities pi come 
out of the extremization procedure as 

Pi = ^{l-{l-q)PEi)Tk, (2) 
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where 

Z, = Y.{l-{l-q)m^^ (3) 

i 

is the partition function, (3 = l/{kT) is a suitable inverse temperature variable, 
and the Ei denote effective energy levels of the microstates i . In the limit g — > 1, 
ordinary statistical mechanics is recovered, and eq. (^) reduces to the ordinary 
Boltzmann factor pi ~ e~^^\ Most formulas of thermodynamics appear to have 
a simple generalization valid for arbitrary q |^]. 

Although this is certainly a beautiful mathematical framework, a natural ques- 
tion is whether this type of statistical mechanics is physically realized. Tsallis' 
original suggestion was that it may have physical applications for equilibrium sys- 
tems with long-range interactions (see |]l9| for a recent update). Very recently, a 
slightly different application was suggested. It was pointed out that the formal- 
ism is of particular physical relevance for nonequilibrium systems with fluctuating 



temperature or fluctuating energy dissipation [|^, 21 1. The basic idea how 

fluctuations can generate nonextensive statistics is easily understood: Consider an 
arbitrary Hamiltonian with energy levels Ei of the microstates. Using the integral 
representation of the gamma function one can easily prove the following formula 

1 

(1 - (1 - q)PoE^)—^ = / e-^^'fimP, (4) 

Jo 

where 

q-l = - (5) 

n 

and 

n 

is the probability density of the (o^' gamma) distribution of order n. It arises 
in a natural way if n independent Gaussian random variables with average are 
squared and summed up. 

Formula @) is valid for arbitrary Ei and thus universal. The left-hand side 
of eq. (^) is just the generalized Boltzmann factor emerging out of nonextensive 
statistical mechanics. The right-hand side is a weighted average over Boltzmann 
factors of ordinary statistical mechanics. If we consider a nonequilibrim system 
with fluctuating temperature, then nonextensive behaviour is a consequence of 
integrating over all possible fluctuating inverse temperatures /3, provided (3 is x^- 
distributed. 

The constant Pq in eq. is the average of the fluctuating /?, 

iP) := / Pf{P)dP = Po (7) 
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and for the relative variance of the fluctuations one obtains 



3 Langevin equation with fluctuating friction 
forces 

We will now construct a concrete dynamics where Tsallis statistics can be rigor- 
ously proved. 

Consider a nonlinear Langevin equation of the form 

ii = - jF{u) + aL{t) (9) 

where F{u) = — ^^(^i) is a nonlinear forcing and L{t) is Gaussian white noise. 
For the turbulence application, we are mainly interested in power-law potentials 
of the form V{u) = C|np". Given some fixed ratio /? := 7/0"^ of the parameters in 
the stochastic differential equation this generates a process with the stationary 
probability density function 

p{u\(3) = iC(3)^ exp {-(3C\u\'"'} . (10) 

^ \2a) 

We now allow the parameters 7 and a in eq. (|9|) to fluctuate as well. Let us assume 
that either 7 or a or both fluctuate in such a way that /5 = 7/(7^ is x^-distributed 
with degree n. For example, we can generate this x'^-distribution by yet another 
set of Langevin equations, say 

Xj = -rXj + Lj{t) (j = l,...,n), (11) 

and by defining 

n 

i=i 

On a long time scale this set of linear Langevin equations clearly generates in- 
dependent Gaussian random variables Xj with average (provided F = const), 
hence P is x^-distributed with density (|6|). 

Assume that the time scale on which f3 fluctuates is very large, so that for a 
given P local equilibrium described by p{u\P) is achieved in eq. (^). In the long- 
term run we have to take into account all fluctuating values of p. This means 
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we will ultimately observe the marginal distribution p{u) := f p{u\f3)f{(3)dp being 
generated by eq. @. This integral can be performed exactly and one obtains 



= -y- z (13) 



where 



and 
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Z,-'=a(c(9-l)g)'° . V"'' .N (") 



\2aJ ^ \g-l 2a 



i + r^ (15) 



2ni 

P = T-r^ /3o. (16) 



an + 1 
2a 

1 + 2a-q 

4 Application to turbulent flows 

Let us now come to our main physical example, fully developed turbulence. Let u 
in eq. (P) represent a local velocity difference in a fully developed turbulent flow 
as measured on a certain scale r. It is well known that the dissipation fluctu- 
ates in turbulent flows, so our model makes sense. The parameter f3 = 7/0"^ of 
the stochastic differential equation is an a priori unknown function of the energy 
dissipation in the flow. Let us consider a model where 

/3 = ^ = e,rA. (17) 

Here is the (fluctuating) energy dissipation rate averaged over and r is a 
typical time scale during which energy is transferred. A is a constant with dimen- 
sion length'^ /time'^, its value is irrelevant for the following. Both and r can 
fluctuate, and we assume that /3 ~ e^r is x^-distributed, thus ending up with the 
dynamical model of the previous section. 

We remind those readers not familiar with turbulence that the averaged energy 
dissipation rate er{x,t) in a volume V of size is defined as 

er(x,t) = ^l e(x + f,t)d^r, (18) 
r-^ Jv 

where 
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is the instantaneous dissipation, v is the kinematic viscosity, and Vi is the velocity 
fluctuation (the deviation from the mean velocity) in i direction [^] . The physical 
idea of considering the generalized Langevin equation (^ with /? ~ e^-T is that 
a test particle in the turbulent flow moves for a while in a certain region with a 
given e^, then moves to another region with another e^, and so on. In this picture, 
we actually associate the fluctuations of f5 with spatial fluctuations in a spatially 
extended system. 

At the smallest scale of the turbulent flow, the Kolmogorov length scale r\ = 
(u^/e)^^^, we can give a theoretical prediction of the entropic index q. Here we 
can write 

P = er^A = {euy/^A = {ur,fA, (20) 

where := [yje)^!'^ denotes the Kolmogorov time and := ?7/r^ = [ye)^!'^ is 
the Kolmogorov velocity (see, e.g., |^). If e fluctuates then all these quantities 
fluctuate as well. Due to the three spatial directions one can actually define 3 
independently fluctuating Kolmogorov velocities u:^, j = 1,2,3 such that = 



Z^i=i(^^)^- The 3 components describe the flow of energy into the 
three different space directions. The simplest model assumption is that these 
Kolmogorov velocities are Gaussian with average 0. This means, we identify Xj = 
VAui) in eq. (|l2|) . Hence at the smallest scale the 3 space dimensions lead to n = 3 
or, using eq. (|T5|), g = | if a = 1. 

Note that our approach allows us to view the fluctuations of e at the Kol- 
mogorov scale in terms of a (hypothetical) ordinary Brownian particle of mass M 
that is subjected to ordinary thermal noise of temperature T. The fluctuating 
velocity V of this constructed Brownian particle coincides with the fluctuating 
vector of Kolmogorov velocities. The Brownian particle just absorbs the tur- 
bulent energy flow at the Kolmogorov scale. It bridges the gap between thermal 
and macroscopic description. Equipartition of energy yields 

lM{{ulf) = \kT. (21) 

Using eq. ( [20| ) we can estimate the mass of our energy- absorbing Brownian particle 
as 

5 Comparison with experiments 

The quality of agreement of the generalized canonical distributions ( [I^ ) with exper- 
imentally measured densities is remarkable. To start with. Fig. 1 shows probability 
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densities of longitudinal velocity differences as measured by Swinney et al. in a 
turbulent Taylor-Couette flow on various scales r |^. They are very well fitted by 
eq. (|l3|) choosing an r-dependent q and a = 2 — q. In in total 152 different 
densities were evaluated, varying the spatial scale r and the Reynolds numbers Re. 
Previous fitting attempts using e.g. stretched exponentials do not yield anything 
of comparable precision. 

Next, Fig. 2 shows the results of measurements in another turbulence experi- 
ment. Here we see a histogram of the acceleration a of a test particle advected by 
the turbulent flow (measurements of Lagrangian turbulence by Bodenschatz et al. 
[23|). The acceleration has been rescaled to variance 1. It can be regarded as a 



velocity difference on the smallest time scale that makes sense in turbulence, the 
Kolmogorov time scale. Hence our consideration in the previous section predicts 



Tsallis statistics with q = 3/2 provided a = 1. Formula (ll3|) yields a very good fit 
of the measured distributions (solid line in Fig. 2). Note that we do not use any 
fitting parameter in this plot, since q = 3/2 is predicted, and variance 1 implies 



/3C = 4. On the contrary, in their paper |23] Bodenschatz et al. fit their data with 
the density function 

p{x) = ci-exp[ ^ I , (23) 

i.e. a somewhat arbitrary modification of a stretched exponential with several 
parameters Cj. It yields a fit of similar quality as our formula but requires much 
more free parameters. Moreover, there is no theoretical reason whatsoever for a 



density of the form (22), whereas nonextensive statistical mechanics, leading to 
eq. dl^), is a well-defined theoretical concept distinguished by an extremization 
principle and generalized Khinchin axioms H]. 

Also scaling exponents can be evaluated using nonextensive methods. The 
scaling exponents Cm of moments of radial velocity differences u = v{x + r) — v{x) 
measured at distance r in fully developed turbulent flows have a long history (see, 
e.g., PI, |2|]). In the inertial range one observes scaling behaviour of the form 



(Inr)-/-. (24) 

Recently it has been suggested to calculate the Cm using nonextensive methods 
[^, 12, 14|. One approach is based on an extension of the multifractal model of 



turbulence |12|-[p.4|, another one on a combination of the approach described here 
with so-called extended self-similarity (ESS) [^^. A typical result of the latter 
approach is shown in Fig. 3 (see Q for more details). The nonextensive model 
clearly reproduces the experimental data of scaling exponents very well. This is, by 
the way, not even surprising, because if the densities (p^) perfectly coincide with 
the experimental densities than the moments must coincide as well. The moments 
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of the generalized canonical distribution ( [T3D can be analytically evaluated, they 
are given by 



- B { "^+1 1 m+l 

{\un = I \urp{n)du = \ — } ^ ' (25) 



>' -° I 2a' o-l 2a 



where B{x,y) denotes the beta function. 

Finally, one may also try to understand the scale-dependence of the entropic 
index q{r), which has been measured very precisely in |Q. Nonextensive statistical 
mechanics does not say anything on the function q{r). For the r-dependence 
one needs additional theoretical input. One possibility was recently suggested in 



[21|. The basic idea is as follows. The Tsallis entropies are non-extensive, i.e. for 



independent subsystems I and II composed to a single system I+II one has |]^, |3|, ^ 

Si + 5f + (1 - q)SlSl' = (26) 

The observation, however, is that it is possible to make the Tsallis entropies quasi- 
additive by choosing different entropic indices at different scales. I.e., given a 
certain q for two small independent subsystems I and II we may choose another q' 
for the larger, composed system I-I-II such that 

S!, + 5f = Sl+''. (27) 

This property is called quasi- additivity [^]. In turbulence, q is close to 1, so that 
a perturbative expansion in q — 1 makes sense. The final result of a calculation in 
pl| is that if higher-order contributions in q — 1 are neglected then quasi-additivity 
implies a power law for l/{q — 1) as a function of r, 

r \ 

(28) 



q-1 \i] 

provided r is large enough. The exponent 6 is non-universal, at least for finite 
Reynolds numbers. Fig. 4 shows that the Swinney data indeed confirm a power 
law. 

To summarize, the nonextensive approach yields a very useful new method to 
describe the most important statistical properties of fully developed turbulence. 
It is in very good agreement with experimental measurements. 
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Figure captions 

Fig. 1 Experimentally measured probability densities in a Taylor-Couette flow 
as measured by Swinney et al. |^ for r = 23, 46, 93, 208, 399, 830?7 and comparison 
with formula (|l|) with q = 1.150, 1.124, 1.106, 1.086, 1.066, 1.055 {a = 2-q). The 
j-th density is shifted by —j units. All densities are rescaled to variance 1. 

Fig. 2 Experimentally measured measured probability density of the acceler- 
ation X := a/(a^)^/^ in Lagrangian turbulence as measured by Bodenschatz et al. 



[23| for Reynolds numbers R\ = 200,690,970 and comparison with formula ([1 



with q = 3/2 and a = 1. 

Fig. 3 Scaling exponents Cm as measured in five different experiments |26]-|3C] 
and comparison with the nonextensive prediction of [^. 

Fig. 4 {q — as a function of r. The data points are the values obtained 
from the the measured densities of Swinney et al. Q, the straight line corresponds 
to a power law as predicted by quasi-additivity. The observed exponent is 5 = 0.30. 
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